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Abstract 



Differential cross sections for elastic scattering of 21.5 MeV positive and 
negative pions by Si, Ca, Ni and Zr have been measured as part of a study 
of the pion-nucleus potential across threshold. The 'anomalous' repulsion 
in the s-wave term was observed, as is the case with pionic atoms. The 
extra repulsion can be accounted for by a chiral-motivated model where the 
pion decay constant is modified in the medium. Unlike in pionic atoms, 
the anomaly cannot be removed by merely introducing an empirical on-shell 
energy dependence. 
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A long-standing 'anomalous' s-wave repulsion in the pion-nucleus interaction at thresh- 
old, as found in phenomenological analyses of strong interaction effects in pionic atoms, 
has recently received considerable attention [1-9] following the suggestion by Weise [1] that 
it could be expected, at least in the isovector channel, to result from a chirally motivated 
approach where the pion decay constant becomes density dependent in the nuclear medium. 
Very recently it was also argued [7,8] that the energy dependence of the chirally expanded 
irN isoscalar and isovector amplitudes bo(E) and b\(E) respectively, for zero- momentum off 
shell pions near threshold, could explain the anomaly. In addition, the empirical on shell 
energy dependence was shown in Ref. [9] to be capable of explaining the anomaly by im- 
posing the minimal substitution requirement [10] of E — > E — V c , where V c is the Coulomb 
potential, on the properly constructed pion optical potential used in a large-scale fit to 100 
pionic-atom data across the periodic table. The availability in recent years of strong inter- 
action level shifts and widths of 'deeply bound' Is and 2p pionic atom states in isotopes of 
Pb and Sn [11-14], obtained from the (d, 3 He) reaction, provides added impetus to the study 
of the s-wave pion-nucleus interaction at threshold. 

In the present Letter we report on extending the study of the s-wave term of the pion- 
nucleus potential by measuring the elastic scattering of very low energy 7r + and 7r~ on several 
nuclei. The purpose of this experiment is to study the behaviour of the pion-nucleus potential 
across threshold into the scattering regime and to examine if the above-mentioned anomaly 
is observed also above threshold. Of particular importance is the question of whether the 
density dependence or the empirical energy dependence, which remove the anomaly in pionic 
atoms, are required also by the scattering data. In the scattering scenario, unlike in the 
atomic case, one can study both charge states of the pion, thus increasing sensitivities to 
isovector effects and to the energy dependence due to the Coulomb interaction. 

It is somewhat surprising to realize that at kinetic energies well below 50 MeV there 
seems to be only one set of high quality data available for both charge states of the pion 
obtained in the same experiment, namely, the data of Wright et al. [15] for 19.5 MeV pions 
on calcium. Focusing attention on the isovector channel, it is desirable to include also nuclei 
with N > Z. Therefore the targets chosen for the present experiment, which was carried out 
at 21.5 MeV pion kinetic energy, were Si, Ca, Ni and Zr, where the last two have an excess 
of neutrons. Natural isotopic mixtures have been used in all cases and have been taken into 
account in the model calculations accordingly. 

The experiment was performed at 21.5 MeV pion energy on the piE3 channel at PSI 
[16], using the high resolution magnetic spectrometer LEPS [17]. Particles were identified 
by time-of-flight relative to the HF signal of the cyclotron and by time-of-flight within the 
spectrometer. Self supporting targets were used in all cases. The beam was monitored 
by four decay-muon telescopes and by a downstream hodoscope [17]. The latter was also 
used to determine the beam composition and profile. The elastic scattering of muons with 
the same momentum as the pions (except for the different energy losses in material such 
as beam line and spectrometer windows and plastic scintillators) was used to calibrate the 
absolute scale of cross sections and to check the overall validity of data reduction. This was 
achieved by comparing the measured angular distributions for elastic scattering of muons 
with predictions made for Coulomb scattering from the known charge distributions of the 
target nuclei. Details on this method of muonic normalization, which has frequently been 
used for low energy pion measurements with the LEPS spectrometer, can be found in Refs. 
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[17,18]. Two measurements of elastic scattering of muons were made for each measurement 
for pions: (i) Muons were recorded in parallel with pions but at a slightly different location 
in the focal plane, due to the different energy losses of muons and pions in the target and 
in the various scintillators, (ii) Muons were recorded in designated muon runs, where the 
spectrometer fields were adjusted to bring the muons to the same location in the focal plane 
as the pions in the proper pion runs. After correcting for effective target thickness and 
detector acceptance we have obtained normalization constants common to all angles and to 
the two types of measurements (i) and (ii) mentioned above, but slightly different for p + 
and p~, due to differences in the muon/pion ratios for the different beams. In this way we 
could confirm the dependence of the acceptance of the spectrometer on the position in the 
focal plane, which was determined separately by scanning the magnetic fields. 

Figure 1 shows as an example comparisons between calculations and measurements for 
the Coulomb scattering of muons by Ni. Open and filled symbols are for the two types 
of muon measurements (i) and (ii) mentioned above, respectively. The reliability of the 
focal-plane position-dependence of the acceptance of LEPS is of major importance as the 
conclusions regarding the pion-nucleus interaction (see below) rely exclusively on the shape of 
the angular distributions. The points at 30° are not plotted because different normalizations 
applied due to different settings of the channel slits. 

Figure 2 shows the experimental results for elastic scattering of pions and predictions 
of best-fit optical potentials. Before discussing the implications of x 2 fits, a comment is in 
order on values of errors. (Full details of the experiment and the results will be published 
elsewhere). The number of pion counts in the elastic scattering peaks was usually greater 
than 3000 and the background in the spectrometer was negligibly small, thus resulting in 
statistical errors of less than 2%. Most of the uncertainties in this experiment come from 
the monitoring of the beam intensity and its composition and from the dependence of the 
acceptance of the spectrometer on the position in the focal plane. The measurements with 
muons provided stringent tests of the latter and careful analysis did not reveal any systematic 
effects. An estimated overall normalization error of 5% was given in [18]. For the purpose of 
X 2 fits we have added 5% in quadrature to the statistical errors for each point individually, 
but in order to check the dependence of the derived parameter values on the errors we have 
repeated the x 2 fits a l so f° r only 3% added in quadrature. Values of fit parameters and in 
particular values of the isovector amplitude b\ which is at the focus of the present work did 
not differ between fits made with these two sets of errors. 

The interaction of pions with target nuclei was described by the Klein-Gordon equation 
with the standard potential due to Ericson and Ericson [19] where double scattering, ab- 
sorption on two nucleons and angle-transformation terms have been included [6,20]. The 
potential is written as 

2pV opt (r)=q(r)+V-a(r)V (1) 
with its s-wave term which is the prime concern in the present work, given by 

q {r) = + j^){b (r)\p n (r) + p p (r)} ± &i[p n (r) - Pp (r)}} 

~ M1 + ^M )ABQpn{r)f)p{r) (2) 
where the ± sign refers to the pion =p charge states, respectively. In these expressions p n and 
p p are the neutron and proton density distributions normalized to the number of neutrons 
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N and number of protons Z, respectively, p, is the pion-nucleus relativistic reduced mass 
and M is the mass of the nucleon. The function b (r) in Eq. (2) is given in terms of the 
local Fermi momentum fcp(r) corresponding to the isoscalar nucleon density distribution: 

b (r)^bo-^(b 2 + 2b 2 1 )k F (r), (3) 

where 60 and b\ are the pion-nucleon isoscalar and isovector effective scattering amplitudes, 
respectively. The parameter Bq represents s-wave absorption on pairs of nucleons. The 
expressions for the p-wave term are given in Ref. [6,20] and will not be given explicitly here, 
except for its linear part, namely, 

a(r) = 4tt(1 + {c [p n (r) + p p (r)} ± Ci[p„(r) - p p (r)}} + quadratic terms. (4) 

Nuclear densities were obtained from charge densities and using for the difference between 
neutron and proton rms radii either the results of relativistic mean field calculations or 
values obtained from antiprotonic atoms, see Ref. [6] for details. Derived values of bi were 
insensitive to assumptions on p n within those limits. 

First attempts at parameter fits to the data using the above potential ran into difficulties 
which could be traced to the two-nucleon absorption in the p-wave term. Subsequent fits to 
the 7r + data and the n~ data separately revealed the need to make the p-wave absorption 
parameter for n~ considerably larger than the corresponding parameter for 7r + , an effect not 
seen in earlier fits to only 7r~ data [21]. This may be expected at 20 MeV since the effects of 
the (3,3) resonance should depend on the energy, which is effectively higher for ir~ than for 
7r + . To avoid introducing isospin dependence into the quadratic p-wave term we dropped 
it altogether and made the parameters c and C\ complex. Recall that only at threshold 
these cannot have imaginary terms. The form of the s-wave part of the potential was kept 
initially as for pionic atoms while its parameters were also varied in the fit process. 

Least-squares fits to the whole 7r + and n~ data together produced reasonably good 
agreement between calculation and experiment, with 60 an d the real parts of cq and c\ 
close to the corresponding free 71N values but with 61 significantly more repulsive than 
the corresponding free ttN value of approximately —0.090 m" 1 . The imaginary parts of 
c and c\ were 0.043±0.009 and 0.45±0.11 m~ 3 , respectively. These must be regarded as 
effective because they now represent all non s-wave absorption processes. The s-wave part 
of the potential which was kept as for pionic atoms showed a factor of two reduction in the 
two-nucleon absorption lmB compared to threshold. This model is denoted below as (a). 

With the value of b\ found to be 'anomalously' repulsive, as in pionic atoms, we then 
applied the two mechanisms which have been found to account for that anomaly in the pionic 
atoms case. The first is that due to Weise [1]: since b\ in free-space is well approximated in 
lowest chiral-expansion order by the Tomozawa- Weinberg expression [22] 

&i = -£^ = -0.08 m; 1 , (5) 

then it can be argued that bi will be modified in the pion-nucleus interaction if the pion 
decay constant f n is modified in the medium. The square of this decay constant is given, 
in leading order, as a linear function of the nuclear density, f* 2 — f 2 — pcr/m 2 with a the 
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pion-nucleon sigma term. This leads to a density-dependent isovector amplitude such that 
bi becomes [3] b\(p) = &i (0) /(l — 2.3p) for <r=50 MeV [23] and with p in units of fm~ 3 . 
This model is denoted by (b). The second mechanism which has been successful in pionic 
atoms is to use the empirical on-shell energy dependence of the scattering amplitude b and 
scattering volume Co within the minimal substitution requirement E — > E — V c of [9,10], 
(model (c)). In additional fits we have included both mechanisms together, (model (d)). 

Table I (upper part) summarizes these fits, where 'mixed' refers to the model where 
the s-wave term of the potential is of the conventional type whereas a linear form is used 
for the p-wave part. In the lower part of the table results are given for fits where the s- 
wave potential, too, was linear; models (a') to (d') correspond to models (a) to (d) above, 
respectively. By replacing the density-quadratic s-wave part by a complex linear part we 
checked the dependence of the conclusions on the model used, particularly regarding the 
in medium value of b±, which is found to be essentially decoupled from the rest of the 
potential. 

The obvious conclusions from the table are that (i)without the density dependence of 
bi, its values differ from the free nN value of —0.090 m^ 1 by m 3-4 standard deviations, 
(ii)with the density dependence included, the values of b\ are in agreement with the free 
7riV value, and (iii)the inclusion of the empirical on-shell energy dependence of b and c 
leads to improved fits to the data. The solid lines in figure 2 show the best fit to the data 
obtained with the empirical energy dependence within the linear model (c'). Applying also 
the density dependence of b\ leads to very similar results. The dashed lines show, as an 
example, results of model (a') for Ca and it is evident that the agreement with the data is 
inferior to that achieved with potential (c'), as is also evident from the table. Note that at 
forward angles there is a limited dependence of the cross sections on the strong interaction 
potential and the agreement between calculation and experiment essentially proves that the 
absolute normalization is right. The information on the optical potential comes from larger 
angles, thus it is the shape of measured angular distributions which is sensitive to the strong 
interaction. If the data of Wright et al. [15] for 19.5 MeV pions scattered from Ca are also 
included in the analysis, then the resulting potential parameters are unchanged but their 
7r~ cross sections appear to be systematically too small at larger angles by typically 10%. 

In conclusion, we have performed precision measurements of elastic scattering of 21.5 
MeV positive and negative pions on targets of Si, Ca, Ni and Zr with the aim of testing 
the 'anomalous' s-wave repulsion observed in pionic atoms. In particular, we focused on the 
question of whether the mechanisms, which were shown recently to be capable of removing 
the anomaly in the atomic case, would be required also in the scattering regime in order to 
reconcile the pion-nucleus interaction with the free pion-nucleon interaction. It is found that 
(i) the in-medium isovector amplitude b\ is too repulsive by 3-4 standard deviations compared 
to the free irN value and (ii) that including the empirical on-shell energy dependence of the 
scattering amplitude b and the scattering volume c improves the fits to the data. However 
and unlike the case with pionic atoms, the energy dependence alone does not remove the 
extra repulsion. Only with the inclusion of the chiral-motivated model where the pion decay 
constant is modified in the medium, it is possible to reconcile the pion-nucleus interaction 
with the free pion-nucleon interaction, thus removing the anomaly in the isovector channel. 
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TABLES 



TABLE I. Results of x 2 fits to the data using potentials discussed in the text. The b\ (p) model 
is given by b\(p) = &i(0)/(l — 2.3p). When 61 is complex (bottom half), the listed values refer to 
its real part. 



potential model 


bo and cq 


bi model 


X 2 for 72 points 


61 (m/) 


mixed (a) 


fixed 


fixed 


273 


-0.125±0.015 


(b) 




h(p) 


273 


-0.098±0.008 


(c) 


b (E) and c (E) 


fixed 


197 


-0.117L0.011 


(d) 




bi( P ) 


215 


-0.087±0.010 


linear (a') 


fixed 


fixed 


309 


-0.138±0.010 


(b') 




bi{p) 


259 


-0.096±0.006 


(C) 


b (E) and c (E) 


fixed 


178 


-0.134L0.010 


(d') 




blip) 


214 


-0.093±0.007 
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FIGURES 
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FIG. 1. Coulomb scattering of muons by Ni. Open symbols: from pion runs. Filled symbols: 
from designated muon runs. Continuous curves are calculated Coulomb scattering for the finite 
size charge distribution. Common normalization constants have been used, separately for all the 
fi + points and all the fi~ points, see text. 
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FIG. 2. Comparisons between experimental and calculated differential cross sections for elastic 
scattering of pions. Solid lines are for the best-fit optical potential (c') of Table I, dashed lines are 
examples for potential (a'), see text. 
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